The stability problem and special solutions for the 5-components 

Maxwell-Bloch equations 
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Abstract 

For the 5-components Maxwell-Bloch system the stability problem for the isolated equilibria is 
completely solved. Using the geometry of the symplectic leaves, a detailed construction of the 
homoclinic orbits is given. Studying the problem of invariant sets for the system, we discover a rich 
family of periodic solutions in explicit form. 

1 Introduction 

After averaging and neglecting non-resonant terms, the unperturbed Maxwell-Bloch dynamics in the 
rotating wave approximation (RWA) is given by 




A*y), 

where A, Y are complex scalar functions, that are denoting the self-consistent electric field and respec- 
tively the polarizability of the laser-matter, Z is a real scalar function, which denotes the difference of 
its occupation numbers. The superscript * stands for the complex conjugate. For more details about 
the history and physical interpretations of this system see [3], [TO], [H]- 

Writing X = x\^- 1x2, Y = yi + 11/2 and Z = z the above system transforms into the 5-components 
Maxwell-Bloch system 

ii = yi 
yi = xiz 

±2 = y2 (1-1) 

7/2 = X2Z 

z = -{xxy\ + X2y2)- 

The Maxwell-Bloch system in the form has the advantage of a rich underlying geometrical structure 
that can be used in the study of its dynamical behavior. 

The system admits a Hamilton-Poisson formulation, where the Poisson tensor is given by 



J{xi,yi,X2,y2,z) 
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and the Hamiltonian function is given by 

H{xi,yi,x2,y2,z) = '^{yl + yl + z^) . 

The system has two additional constants of motion, namely the Casimir of the Poisson structure J, 
given by 

C[xi,yi,X2,y2,z) = -{x\ + xj) + z 
and a constant of motion derived from a bi-Hamiltonian structure of the system (see |10] ) given 

by 

I{xi,yi,X2,y2,z) = X2J/1 - xiy2. 

A commuting property of the constants of motion H and / holds, i.e. {H,I} — 0, where {•, •} is the 
Poisson bracket associated to the Poisson tensor J. 

2 Stability of equilibria 

By a direct computation we obtain three families of equilibria for the system 

£1 = {(0,0,0,0,M)| Af e M*}; £2 = {(A/, 0, A^, 0, 0)| A/,iV e R, + ^ 0}; £3 - {(0,0,0,0,0)}. 

It is a well known fact that the dynamics of a Hamilton-Poisson system is foliated by the symplectic 
leaves associated to the Poisson structure. In our case the regular symplectic leaves are given by the 
connected components corresponding to pre-images of regular values of the Casimir function C. We 
denote by Oc — C^^(c),c e R the regular symplectic leaves of the Poisson structure J. 

The restriction of the dynamics (|l.ip to a regular leaf Oc becomes a completely integrable Hamilto- 
nian system 

(0„c^o.,if|oJ, (2.1) 

where the second commuting constant of motion is I\oa- We will study the stability problem of equilibria 
on regular leaves Oc analogously to the approach used in 2 . 

The equilibria of the Hamiltonian system (12.11) can be divided in two types: 

Xq := {{xi,yi,X2,y2,z) e Oc \ d{H\oJ (xi,yi,x2,y2, z) = 0, d(/|oJ {xi,yi,X2,y2, z) 0}; 
Xi := {{xi,yi,X2,y2,z) e Oc | d{H\oJ {xi,yi,X2,y2, z) = 0, d(/|oJ {xi,yi,X2,y2, z) ^ 0}. 

Proposition 2.1. On a reyular symplectic leaf Oc we have the following characterization for the equi- 
libria: 

3Co-Ocn(£iU£3); 3Ci = Ocn£2. 

Proof. Because (|2.ip is a Hamiltonian system on a symplectic manifold the condition d {H\o^) (e) = 
is verified for any equilibrium point e G Oc- 
Let 62 e Oc n £2. Then 

Tc.Oc = {C = (Wl, V2, W3, «4, «5) e R'^l < V: VC(e2) 0} = {w G M^j ViM + + V5^ 0}. 

We also have d/(e2) = Ndyi - Afdy2- Taking, for example, v = {~N,N,M,-M,Q) e Te^Oc we have 
d(/|oJ (e2)(w) = A'P ^ 0, which proves that d(J|oJ (es) 0. 

For the equilibria e in £i U £3 the condition d (/|oc) (g) = is trivially verified. □ 

The commutativity of the constants of motion H\o^ and /|oe with respect to the symplectic form 
ujo^ implies that at an equilibrium point e G Oc we have 

[DXff|„Je),DX,|„Je)] =0, 

where IDX^^^ (e) and DXi^^ (e) are the derivatives of the vector fields Xj^^^ and Xj^^ at the equilib- 
rium e and consequently DXfj^^ (e), T)Xj^^ (e) are infinitesimally symplectic relative to the symplectic 
form ujo^{e) on the vector space TeOc- 
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Definition 2.1. An equilibrium point e € OCq is called non- degenerate i/DX^|g (e) and DXj^^ (e) 
generate a Cartan subalgehra of the Lie algebra of infinitesimal linear transformations of the symplectic 
vector space (TgOc, wo^(e)). A Cartan subalgebra of the Lie algebra sp(4,]R) is a two dimensional 
commutative sub-algebra which contains an element whose eigenvalues are all distinct. 

It follows that for a non-degenerate equilibrium belonging to 3Co the matrices DX//|g^(e) and 
DX/|^ (e) can be simultaneously conjugated to one of the following four Cartan sub-algebras 



Type 1: 
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(2.2) 



where A,B£M. (see, e.g., 0], Theorems 1.3 and 1.4). 

Equilibria of type 1 are called center- center with the corresponding eigenvalues for the linearized system: 
iA, -lA, iB, -iB. 

Equilibria of type 2 are called center-saddle with the corresponding eigenvalues for the linearized system: 

A,~A,iB,-iB. 

Equilibria of type 3 are called saddle-saddle with the corresponding eigenvalues for the linearized system: 
A,-A,B,~B. 

Equilibria of type 4 are called focus-focus with the corresponding eigenvalues for the linearized system: 
A-\-iB,A- iB, -A + iB, -A - iB. 

Theorem 2.2. We have the following stability behavior for the equilibria in Oc H £i.' 

(i) The equilibrium point Oc H £i ~ {(G, 0, G, 0, c)} for c > is a non- degenerate equilibrium of type 
focus-focus and consequently unstable. 

(ii) The equilibrium point Oc H £i = {(G, 0, G, 0, c)} for c < is a non- degenerate equilibrium of type 
center-center and consequently stable. 

Proof, (i) For the linearized systems at the equilibrium (G, G, 0, G, c) we have: 

DXh|„JO,G,0,G,c) = 

and its characteristic polynomial has the non-distinct eigenvalues y^, -y/c, — y^, —\/c and respectively 

DX,|„jO,G,0,G,c): 



1 G G 

c G G 

G 1 

c G 
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-1 

-1 
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1 







and its characteristic polynomial has the non-distinct eigenvalues i, i, —i, —i. 

To decide the type of stability we need to determine the non-degeneracy of the equilibrium (G, G, G, 0, c), 
i.e. we have to find a linear combination TiXjj^^ (G,0,0, G,c) -\- aT>Xj\^ (0,0,0, G,c), where a is a 
non-zero real number, that has distinct eigenvalues. The characteristic polynomial of such a linear 
combination is given by 

t^ + {2a^ -2c)t' -^{a" +cf. 
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After the substitution = s we obtain the quadratic polynomial 

+ {1o? - 2c)s + (a^ +c)2, 

which has the discriminant A ~ —\^co? < and therefore has two distinct complex roots. It follows 
that the characteristic polynomial t^ + (2a^ — 2c)t^ + (a^ + c)^ has four distinct complex eigenvalues of the 
form A + iB,A-iB, -A + iB,-A - iB with A,B€R*. Consequently, the equilibrium (0,0,0,0, c) for 
c > is a non-degenerate equilibrium of focus-focus type for the dynamics (|2.ip and therefore unstable 
for this dynamics. 

Similar computations lead to the proof of (ii) . □ 

Although the equilibrium (0, 0, 0, 0, 0) G Oq belongs to 3Cq, it is a degenerate equilibrium in the sense 
of Definition Indeed, any linear combination aBXnia^ (0, 0, 0, 0, 0) + /3DX/|„^ (0, 0, 0, 0, 0) has the 
characteristic polynomial {t^ + which has non-distinct eigenvalues. Its stability property can be 

established using an algebraic method (see [T], [5], [S], [7]). More precisely, the system of algebraic 
equations 

i?(xi, 2/1, X2, 2/2, z) = H{0, 0, 0, 0, 0), /(xi, 2/1, 3^2, 2/2, z) - /(0, 0, 0, 0, 0), C(a;i, yi, X2, y2,z)^ C{0, 0, 0, 0, 0) 
has as unique solution the equilibrium (0,0,0,0,0), leading to the following stability result. 
Theorem 2.3. The equilibrium (0,0,0,0,0) is degenerate and stable with respect to the dynamics (11.11) . 



3 Homoclinic orbits 

In this section we will give an explicit form of the homoclinic orbits for the unstable equilibria of 
focus-focus type. This type of equilibria belong to symplectic orbits Oc with c > 0. 

In order to compute the homoclinic orbits, we introduce a local system of coordinates around the 
equilibrium point Cc = (0, 0, 0, 0, c) G Oc. The local system of coordinates is given by 



Xl 
X2 

yi 

z = c 



ri cos 9 
ri sin 9 

yi 

y2 

1„2 
2^1- 



Freezing the parameter c we obtain the local system of coordinates on the symplectic orbit Oc around 
the equilibrium point Cc- 



Oc\{ec}, (ri, 61,2/1, 2/2) ^ 



Xl 


= Ti COS 


X2 


— ri sin 


yi 


= yi 


y2 


= 2/2- 



As we have excluded the equilibrium point Cc we can work under the assumption that ri ^ 0. The 
advantage of using polar coordinates in the study of bi-focal homoclinic orbits in four dimensions can 
be ascertained in [S]. By a straightforward computation we obtain that the reduced system on the 
symplectic leaf Oc is given by 

ri = 2/1 cos 9 + 1/2 sin 9 
■ _ 2/2 cos 6> - 2/1 sin 9 

(3-1) 

2/1 = ri cos 9 (c - 
2/2 = ri sin 9 (c - 



I 



Using a continuity argument and the fact that I{xi, 2/1, X2, 2/2, z) = X2yi — 2:12/2 is a constant of motion 
we obtain that if there exists a homoclinic it should belong to the connected component of level set 
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I ^(/(cc)) = / ^(0) that contains Cc- If a curve c(t) = {ri{t),9{t),yi{t),y2{t)) is a homoclinic, then it 
has to be a solution for the system p.ip and to satisfy the following equation for all t: 

ri{t){yi{t)sm9{t) - cos0{t)) = 0. 

This implies that 9{t) — and thus 9{t) — do constant for all t. By differentiation and substitution we 
obtain the following second order equation 

h = r-i ^c- ^rj^ . 

Making the change of variable ri = 2-\/c fi and the time re-parametrization ^/c t = t we obtain the 
equation 

fi{i) = fi{i)-2fl{i). 

It is well known that this second order differential equation has as solutions ±cn(<, 1) — ±sech(i). 
Consequently, we obtain ri{t) — ±2-^0 sech(-yc<). 

Substituting ri{t) in the expression of the local parametrization $c, and for z in the expression of 
local parametrization $ and integrating for j/i and 7/2 in (|3.1I) we obtain the homoclinic solutions 

xi{t) — ±2-^0 sech(-yct) cos 6*0 
X2{t) — ±2-^0 sech(Y^t) sin^o 
< yi{t) — t2c sech(Y^t)tanh(-yct) cos 6*0 
y2{t) = =f2c sech(-y/ct)tanh(-yct) sin^o 
= c(l - 2sech^(yci)). 

The above homoclinic orbits, using different parametrization and arguments, have been discussed in [5], 

m- 



4 Invariant sets and periodic orbits 

We will look for invariant sets of the system using the technique presented in |3|. We have the 

following vectorial conserved quantity F : R^,F(p) = {H (p) , I (p) , C (p)) . In [3], Theorem 2.3, it 

has been proved that the set M^^) = G rank VF(p) = 2} is invariant under the dynamics of the 
system. By direct computation we obtain that M^^j ~ Mi U M2, where 

Ml := |(^a;i,yi,X2,-^^^,-%) \ X2 ^ 

I V X2 X^J 

M2 1(^X1, 0,0, y2,-^) \xi^O 

The union Mi U M2, which is a connected set in R^, is invariant under the dynamics but neither 

the set Ml nor the set M2 are invariant under this dynamics. The vector field corresponding to (jl.ip is 
tangent to the sub-manifold Mi and the restricted dynamics on Mi is given by 

xi = yi 

(4.1) 

We notice that the above dynamical system has two conserved quantities, /i , /2 : Mi — > M, fi{xi,yi, X2) = 
x\ + x\ and f2ixi,yi,X2) = f^. Using these conserved quantities and choosing an initial condition 
Xi,yi, X2 with X2 ^ and j/J ^ we can exphcitly solve the system (|4.ip : 



xi{t) = x° sin + a^? cos 

W = Bin cos gi) 

X2(t) = -x1 Sin (^t) + cos . 
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Notice that if = '^^ obtain as constant solutions the equihbrium points from £2 C Mi. The above 
solution is defined on time intervals {tk,tk+i), where tk — -^"i? + ^^'^ with fc G Z and G [0, 27r) is 



the unique real number such that sin 



and cos -d — 



= . For time values tk the 



above solution exits the set Mi and punctures the set M2 , thus making the union Mi U M2 an invariant 
set. Although the solution starting from Mi is not complete, we can construct a complete periodic 
solution for the initial system (jl.ip given by 

xi{t) = x§ sin + xl cos 

,i(t) = -||(.?;in(||i)-xOcos(||i) 
X2it) = -xl sin (||<j + x°2 cos 
y2W--S(x°sin(f|i)+x?cos (gi) 
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